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ABSTRACT 

There is growing evidence that relativistic jets in active galactic nuclei undergo extended (parsec-scale) ac- 
celeration. We argue that, contrary to some suggestions in the literature, this acceleration cannot be purely 
hydrodynamic. Using exact semianalytic solutions of the relativistic MHD equations, we demonstrate that 
the parsec-scale acceleration to relativistic speeds inferred in sources like the radio galaxy NGC 625 1 and the 
quasar 3C 345 can be attributed to magnetic driving. Additional observational implications of this model will 
be explored in future papers in this series. 

Subject headings: galaxies: active — galaxies: individual (NGC 6251) — galaxies: jets — ISM: jets and 
outflows — MHD — quasars: individual (3C 345) 



1. INTRODUCTION 

Magnetic acceleration and collimation has long been 
thought to be the underlying mechanism responsible for the 
similar manifestations of cosmic jets in such diverse systems 
as y oung stellar objects and active galactic nuclei ( AGNs) 
(e.g.. lKonig ]|[l986: Pringle1 [l993l ISpruitl[l996t IUviol l2000). 
Although some AGN jets have long been known to ex- 
hibit apparent superluminal motions, with inferred (termi- 
nal) bulk Lorentz factors in the blazar class of sources of 
7oo 5; 10 (but exceeding 40 in some cases; e.g., Jorstad et al. 
2001), most models to date have concentrated on the non- 
relativistic regime. However, two recent discoveries — 
the detection of apparent superluminal motions in certain 
Galact ic black-h ole binaries (the so-called microquasars; e.g., 
Mirabel & Rodriguez! 1 1 9991) . from which mildly relativistic 
bulk velocities have been deduced, and the inferred associa- 
tion of gamma-ray bursts (GRBs) with ultrarelativistic (700 > 
10 2 ), highly collimated outflows (e.g.. lPirar3ll999h — have 
highlighted the strong similar ities among the various t ypes 
of relativistic jet sources (e.g., Ghisellini & Celotti 2002) and 
have refocused attention on the question of their origin. Al- 
though the interpretation of relativistic outflows and the gen- 
eralization of magnetohydrodynamics (MHD) to the relativis- 
tic regime present several distinct challenges, the prevailing 
view has been that magnetic driving is the common u nderly- 
ing mechanism also in this case (e.g., Blandford 2002). 

However, this interpretation of relativistic jets is by no 
means universal. One example is p rovided by the radio galaxy 
NGC 6251, in which ISudou et alJ feOOOl) inferred a bulk ac- 
celeration from V « 0. 1 3 c to V ss 0.42 c on sub-parsec scales . 
This behavior was attributed bv lMelia. Liu. & FatuzzoM2002l) 
to a thermal acceleration of a proton-electron plasma that 
is heated to a temperature T ps 10 12 K in a region of radius 
r < 0.03 pc. We note, however, that a thermally driven, purely 
hydrodynamic flow typically undergoes the bulk of its ac- 
celeration over a distance that is of the order of the size of 
the mass distribution that initially confines it by its gravity, 
which in this case is much smaller than the radius of the ap- 



parent acceleration zone ( see §|3. Centrifugal driving (e.g., 
Blandford & Pavne 1982) — the commonly invoked hydro- 
magnetic acceleration mechanism for nonrelativistic jets — 
typically also acts fairly rapidly and thus would similarly fail 
to account for the large-scale acceleration inferred in NGC 
6251. A possi ble resolution of this puzzle is provided by 
the finding of Li, Chiueh, & Begelman ( 1992), who were the 
fir st to generalize the "cold" r adially self-similar MHD flows 
of Blandford & Pavne ( 1982) to the relativistic regime (see 
also IContopoulosl l 19941) . that their solutions contain an ex- 
tended magnetic pressure-gradient acceleration region beyond 
the classical fast-magnetosonic point (a singular point of the 
Bernoulli equation), a behavior that they ascribed to the ac- 
tion of a "magnetic n ozzle." It was subsequently shown by 
Vlahakis et al. (2000) that a similar mechanism operates also 
in nonrelativistic flows, but the effect is probably easier to dis- 
cern observationally in relativistic jets. 

Vlahakis & Konigl (2003a, hereafter VK) carried out a 
further generalization by deriving "hot," radially self-similar, 
relativistic MHD solutions for trans-Alfvenic flows. 1 They 
showed that the magnetic field always guides and collimates 
the flow, but that, if the specific enthalpy £c 2 is initially 
(subscript 2) ^> c 2 , then an extended thermal acceleration re- 
gion can develop, within which the flow is accelerated from 
ji ps 1 to 7 w If the total energy-to-mass flux ratio /ic 2 
is £,-c 2 , corresponding to a Poynting flux-dominated out- 
flow, then the bulk of the acceleration is magnetic and takes 
place downstream from this point. VK demonstrated that 
the flow continues to be accelerated all the way up to the 
modified fast-magnetosonic surface, which is the locus of the 
fast-magnetosonic singular points of the combined Bernoulli 

1 The trans-Alfvenic solutions correspond to a dominant poloidal magnetic 
field at the base of the flow. Vlahakis & KonieJ 1 2003b) derived analogous 
solutions for super-Alfvenic jets, for which the magnetic field at the base 
of the flow is predominantly azimuthal. The latter configuration m ay be ex- 
pected to apply to inherently nonsteady outflows (e.g., ContoDoulos 1995). In 
this paper we adopt the poloidal field configuration as the most appropriate 
modeling framework for AGN jets, as has also been done by other workers 
fe.g.. lLivio. Pringle. & Kind20im 
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and transfield equations and represents the true "causality sur- 
face" (or "event horizon") for the propagation of fast waves. 
VK showed that this singular surface can lie well beyond the 
classical fast-magnetosonic surface and argued that this is the 
essence of the "magnetic nozzle" effect. 

This is the first in a series of papers in which we apply 
the VK formalism to the interpretation of relativistic jets in 
AGNs. 2 Our aim is to model a variety of observational find- 
ings and attempt to construct basic diagnostic tools for the 
study of such jets. In this paper we focus on the extended- 
acceleration signature of magnetically driven jets and use our 
solutions to model the parsec-scale accelerations already in- 
dicated in a number of relativistic jet sources. We first present 
arguments for why the extended acceleration is unlikely to 
have a thermal origin (§ |2j. We then consider magnetic jet 
models: after a brief review of the solution methodology 
(§121, we demonstrate (§0} that magnetic driving can account 
for Sudou et al.'s observations of NGC 625 1 as well for the 
parsec-scale acceleration to > 10 inferr ed in superluminal 
blazar jets like 3C 345 ( Unw inet al.ll 1 99% . Our conclusions 
are given in §|5] 

2. CAN THE JETS BE DRIVEN THERMALLY? 

The sub-p arsec-scale k inematics of the NGC 625 1 jet was 
deduced by ISudou et al.l ( 120001) after they discovered the 
counterjet in this powerful radio galaxy. By applying a rel- 
ativistic beaming model to the measured variation in the jet— 
counterjet intensity ratio, they inferred that the outflow is ac- 
celerated from ~ 0.13 c at r w 0.53 pc to ~ 0.42 c at rsj 
1 . pc. 3 We now arg ue that, contrary to the suggestion made 
in lMelia et al.l ll2002l) . the cause of this acceleration cannot be 
thermal pressure driving. 

Since the inferred speeds are not too close to c, it is suf- 
ficient to use nonrelativistic hydrodynamics. In the absence 
of magnetic forces and assuming spherical symmetry for sim- 
plicity, the flow can be described as a Parker ( 1958) wind, in 
which the radial velocity V is given as a function of the radius 
r by 

r dV _ 2C 2 -GM/r 

v~d7~ v 2 -c 2 ' ( } 

Here M is the mass of the central black hole, G is the gravi- 

1 /2 

tational constant, and C s = (iTIc^T /nip] is the sound speed 
of a fully ionized hydrogen gas (with T, m p and k B being the 
adiabatic index, proton mass, and Boltzmann's constant, re- 
spectively). In the best-fit model of Meli a et all ( 120 02) (which 
is consistent wit h the limits set by the radio observations of 
iJones et alii 19861) . the base temperature (at r, = 0.026 pc) is 
Tj = 10 12 K. For these values of r a nd T, and with M as 
6 x 10 s M JFerrarese & Fordl fl999l) . we find that 2C 2 is 
^> GMj r and hence that the flow at r, would already be super- 
sonic with V > C s s» 0.55 c. By the time such a flow reached 
the scales mapped by Sudou et al. (2000), its velocity would 
be significantly larger — and its acceleration substantially 
smaller — than what has been inferred from the observations. 

2 We have previously concentrated on applications to GRBs; see VK, 
where analogies among GRBs, AGNs, and microquas ars are discussed, as 
well a s Vlahakis & Konigl (2001, 2003b) and Vlahakis. Peng. & Konigl 
(2003). 

3 Using a more accurate determination of the distance to NGC 6251, 
IMelia et al] 1200211 took the radial range of the acceleration region to be 
0.30 — 0.57 pc. Since the differences from the values adopted by Su dou et alj 
1 2000 ) have little impact on our arguments, we continue to use the latter in 
our discussion. 



Extended acceleration has also been indicated in super- 
luminal blazar jets. For example, VLBI images of the 
quasar 3C 345 have shown that the jet co mponent speeds in- 
crease with separation fro m the core (e.g.. lZensus et alll995l 
Loban ov"& Zensi^l LL9_9_9j). In par ticular, in the case of the 
C7 component, Unwin et al. ( 1997) combined a VLBI proper- 
motion measurement with an inference of the Doppler factor 
from a synchrotron self-Compton calculation to deduce an ac- 
celeration from 7 ^ 5 to 7 > 10 over a (deprojected) distance 
range (measured from the core) of ~ 3 - 20 pc. 

Although blazar jets are also sometimes modeled 
in terms of a purely hyd rodynamic acceleration (e.g. 
Georganop oulos & Marschen I1998I) . one can show quite 
generally that in this case, too, thermal driving alone can- 
not account for the observations. A variety of arguments 
(summarized, e.g., in Ghisellini & Celotti 2002) indicate that 
protons are the dynamically dominant component in many 
AGN jets. Energy conservation in such outflows implies 
7oo/7; ~ 1 + 2[T/(T- l)](k B Ti/m p c 2 ). Even if the initial 
temperature were as high as ~ 10 12 K, the terminal Lorentz 
factor would still be < 2 — much smaller than the values 
typically inferred from the observed superluminal motions. 
An additional argument can be made on the basis of the 
observed acceleration rate, which in a purely hydrodynamic 
model is determined by the radial dependence of the external 
pressure that provides lateral confinement of the jet. In a 
simplified picture of a one-dimensional pressure distribution 
that scales with distance z from the origin as z~ a , the Lorentz 
factor of an adiabatic, supersonic jet is predicted to increase 
as z"/ 4 (e.g., Blandford & Ree slll974]). Applying this picture 
to the 3C 345 data reported in Un win et all ([1997), one infers 
a ~ 1.5 on parsec scales. It is, however, unclear how a 
pressure distribution of this type could arise in a natural way 
in an AGN. 

The most likely alternative in both of these cases is mag- 
netic acceleration, which we consider in the remainder of this 
paper. 

3. SEMIANALYTIC SOLUTIONS OF THE "HOT" RELATIVISTIC 
MHD EQUATIONS 

VK considered the full set of special-relativistic MHD 
equations, allowing both the bulk and the random speeds to 
be relativistic. They assumed ideal MHD, axisymmetry, no 
explicit time dependence, and a polytropic equation of state 
(with the pressure P scaling with the rest-mass density po as 
P oc pQ, where the adiabatic index is taken to be 5/3 or 4/3 de- 
pending, respectively, on whether the pressure is dominated 
by "cold" protons and electrons or by "hot" electron-positron 
pairs and radiation). Under these assumptions, the MHD 
equations can be partially integrated to yield several field-line 
constants (with the field line being identified by the poloidal 
magnetic flux function A): the total specific angular momen- 
tum L(A), the field angular velocity fl(A), the magnetization 
parameter om(A) (with the mass-to-magnetic flux ratio given 
by Ail 2 /a M c 3 ), the adiabat Q(A) = P/pf ) , and the total energy- 
to-mass flux ratio p(A)c 2 = ^^-(c/AirXEB^/jpoVp) (where 
E is the electric field amplitude, is the azimuthal field com- 
ponent, and V p is the poloidal velocity component). For the 
adopted equation of state, the specific enthalpy is given by 

£c 2 = c 2 +[r/(r-i)]CP/ A) ). 

The transfield force-balance equation is integrated under 
the most general ansatz for radial self-similarity [in spherical 
coordinates (r , 9 , </>))], in which the shape r(A , 9) of a poloidal 
field line is given as a product of a function of A times a func- 
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tion of 9: r= Fi(A)F0) (see IVlahakis & Tsinganoil998l) . 
Such a separation of variables can be effected if !F\(A) oaA l ' F , 
L(A) oc A l ' F , fi(A) (xA~ l ' F , Q(A) (x A~ (F ~ 2)(r ~ l \ and if n(A) = 
const and <Jm(A) = const. To obtain a solution, it is neces- 
sary to specify seven boundary conditions as well as the val- 
ues of r and F. The parameter F controls the distribution of 
the poloidal current /: 2I/c= (SB^ = A l ~ x l F T(0) [using cylin- 
drical coordinates (€3,0,z)]. Close to the origin the field is 
force-free, with T(9) ~ const, which implies O3B0 oc A l ~ l l F . 
For F > 1, the current |/| is an increasing function of A, 
corresponding to the current-carrying regime (for which the 
poloidal current density is antiparallel to the field). A solution 
with F > 1 should provide a good representation of the condi- 
tions near the axis of a highly collimated flow. Conversely, a 
solution with F < 1 corresponds to the return-current regime 
(for which the poloidal current density is parallel to the field) 
and may be most suitable at larger cylindrical distances. VK 
showed that, for F > 1, the modified fast surface is at infin- 
ity and an initially Poynting-dominated flow (p ^> £,) attains 
a rough equipartition between the kinetic and Poynting en- 
ergy fluxes at large distances from the origin. In this case the 
Lorentz force efficiently collimates the flow, which reaches 
cylindrical asymptotics. In contrast, for F < 1, even though 
the acceleration is more efficient (so more of the Poynting flux 
is converted into kinetic energy), the collimation is weaker 
and the flow only reaches conical asymptotics. 

The illustrative solutions derived in VK demonstrated that 
centrifugal driving plays a limited role in the acceleration 
of relativistic flows — it is only important initially (for as 
long as Vc/, > V p ). As already noted in § [2 thermal accel- 
eration can dominate over a more extended zone if the ini- 
tial enthalpy is relativistic (£,■ ^ 1). Magnetic acceleration, 
however, takes over after £ drops to ~ 1 and can remain im- 
portant well beyond the classical fast-magnetosonic surface 
(where jV p « [(B z - E 2 )/(4irp O] l/2 )- If the Poynting flux 
is initially smaller than the enthalpy flux (corresponding to 
[i w £,), the acceleration is predominantly thermal (governed 
by £7 ss fi) and, for a trans-Alfvenic solution, the magnetic 
field only acts to guide the flow. The collimation in this case 
is w eak, and the streaml ines are asymptotically conical (see 
also lVlahakis at al.l2003l) . 

4. MAGNETIC OUTFLOW MODELS FOR JETS IN RADIO GALAXIES 
AND QUASARS 

4.1. The Accelerating Jets in the Radio Galaxy NGC 6251 

Figure shows a solution of the steady, axisymmetric, 
ideal-MHD equations that describes an outflow from a disk 
around a supermassive black hole. As in VK, the solution 
was constructed using the r self-similarity formalism but was 
restricted to extend over a finite radial range. The dimensional 
parameter values were determined on the basis of the inferred 
magnitudes of the NGC 625 1 black-hole mass and mass out- 
flow rate. 4 Given that the slow-magnetosonic singular surface 
arises from the interplay between gravitational and thermal 
forces and that the relativistic r self-similar model does not 
incorporate gravity, the solution only covers the super-slow 
regime of the flow. This regime, however, is the most perti- 
nent one and — even for the high value of 7) that we adopt 

4 The values of the model parameters and boundary conditions that define 
the displayed solution are, in the notation of VK: T = 5/3, F = 1.0001, Zc = 
0.7 pc, x\ = 0.87, cr M = 3, £ A = 1.005, n = 7.8, 6» A = 50°, and Bo®^ = 
1.1 X 10 16 cgs. 



(following iMelia et alJl200l — contains the V „ > 0.13 c ve- 
locity range measured in the Sudou et al. (2000) observations. 

Figure^ depicts the field-line shape. The innermost field 
line originates from the vicinity of the black hole, at a (cylin- 
drical) distance 03i n of a few Schwarzschild radii from the cen- 
ter (for a black hole of mass M, the Schwarzschild radius is 
~ 10~ 13 (M/Mq) pc), whereas the outermost field line orig- 
inates at 03 out = 10 3 03; n . For this choice of the 03 out /03i n ratio 
and the inferred mass of the NGC 625 1 black hole, the mass- 
loss rate in the wind is ~ 9 x 10 21 g s" 1 ( see Fig. \l]p), close to 
the best-fit estimate of lMelia et alJ (120021) . 

Figure shows the various force densities in the poloidal 
direction along the outermost field line as functions of 03/03 a 
(where 03 a = (L/fifl) 1 ^ 2 is the Alfven lever arm; note that, in 
the r self-similar model, 05/03 a is solely a function of the polar 
angle 9). Although the thermal pressure gradient is the domi- 
nant force density very close to the origin, the magnetic force 
rapidly takes over. (The centrifugal force, which could in 
principle also contribute near the origin, is much smaller than 
the pressure gradient force in this case.) FigureQJl shows that 
the Lorentz factor (solid line) increases monotonically with 
distance. Also shown is the product £7 (dashed line), which 
demonstrates (in accord with Fig. ^;) that £ differs from 1 
(signaling that thermal effects contribute to the acceleration) 
only in the immediate vicinity of the origin (03 < 0.35 03a). 
It is seen that, even for very high initial temperatures (Fig. 
[if), the thermal contribution to the acceleration is negligible, 
confirming the conclusion of § [2] The upper curve in Fig- 
ureQJl indicates that, for 03 > 0.35 03a, the Lorentz factor in- 
creases due to a decreasing Poynting-to-mass flux ratio: in 
this range the Poynting flux is converted into matter kinetic- 
energy flux. Asymptotically the Lorentz factor is ~ /i/2, or, 
equivalently, the flow reaches a rough equipartition between 
Poynting and kinetic energy fluxes. Although the final value 
of 7 in the jets in NGC 6251 is not known, our particular 
choice of value for /1 implies that 700 w 4. This estimate is 
consisten t with the values typic ally inferred for radio-galaxy 
jets (e.g., Gio vannini et all2001l) . 

Figure^ shows the various components of the flow ve- 
locity. Although it is not crucial for a qualitative analysis, 
we note that the solution presented h ere act ually reproduces 
the velocities inferred by Sud ou et afl l l2000i) . Asymptotically 
Vm <C V z , as expected in solutions that exhibit cylindrical col- 
limation (see §0. 

The components of the magnetic field along the outermost 
field line are shown in Figure [Q*. The poloidal field initially 
exceeds the azimuthal component (with the two becoming 
comparable at 03 = 03a), but at larger distances the azimuthal 
component dominates. The fields scale roughly as B^ oc 05" 1 
and B- oc 03" 2 . 

4.2. The Accelerating Jet in the Quasar 3C 345 

We propose that the parsec-scale acceleration inferred by 
lUnwin et alJ (11997 ) in component C7 of the 3C 345 jet is most 
plausibly interpreted in terms of magnetic driving, and we 
present in Figure |2] an r self-similar MHD solution describ- 
ing a proton-electron jet that supports this claim. 5 In particu- 
lar, it is seen from Figure |2jl that the Lorentz factor increases 
from 7 ps 5 at 03/03 A ~ 5.55 to 7 « 10 at 03/03 A ~ 13.31. 

5 The parameters/boundary conditions that determine this solution are, in 
the notation of VK: r = 5/3, F = 0.99, Zc = 0, x\ = 0.987, cr M = 8.66, £ A = 
1.0007, ju = 75, 9 A = 25°, and B tBo" F = 2.1 X 10 17 cgs. 
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tD/G3 A G5/G5 A 

FIG. 1 . — r self-similar solution describing the jets in NGC 6251. (a) Poloidal field-line shape on a logarithmic scale, (b) Mass loss rate as a function of 
G5out/®iii! tne rat '° °^ me ou t ermos t and innermost disk radii. The remaining panels show the force densities in the poloidal direction (c) and various other 
quantities [(d) — (g); see text for details] as functions of OS/OSa (which, in turn, is a function of the polar angle 9) along the outermost field line. Here G5a is 
the Alfven lever arm, which equals Ha. in = 9.8 X 10 -4 pc and GJAout = 10 3 ^A.in = 0.98 pc on the innermost and outermost field lines, respectively. Along the 
innermost field line, the quantities plotted in panels (d), (e), and (f) remain the same functions of fTi/fDA, whereas the quantities shown in panels (c) and (g) are 
G3A,out/ti3A,in times larger. 



Using the value of the Alfven lever arm on the outer field 
line (B5A,out = 4.1 x 10~ 2 pc), we find that the cylindrical dis- 
tance changes from G3 as 0.23 pc to 03 « 0.546 pc; on the basis 
of Figure [2^, these cylindrical radii correspond to linear dis- 
tances from the origin of ~ 3 pc and ~ 20 pc, respectively, in 
close correspondence with the observed values. 

The quasar jet solution shown in Figure [2] is characterized 
by a significantly higher value of the total energy-to-mass flux 
ratio /ic 2 than the radio-galaxy solution depicted in Figure 
^ The terminal Lorentz factor, which again corresponds to 
a rough equipartition between the asymptotic Poynting and 
kinetic energy fluxes (700 w /i/2) is correspondingly higher. 
For the adopted fiducial parameters, the C7 component is pre- 
dicted to continue accelerating up to w 35. Interestingly, 
Lorentz factors of this order have been inferred in the more 
dist ant components (in par ticular, C3 and C5) of the 3C 345 
jet tLobanov & Zensusll9 99). 

In comparison with the radio-galaxy solution of § 14.11 the 
quasar solution presented in this subsection corresponds to 
a more massive outflow (with a mass-loss rate ~ 10 23 g s -1 
for G3 0llt /G3j n = 150; see Fig. |2ji) and to a stronger magnetic 
field (cf. Figs. [Q; and|2^). We note, however, that the den- 



sity and magnetic field strength are not uniquely determined 
from the kinematic data: exactly the same flow speeds and 
field-line shape are obtained if the density, particle pressure, 
and squared amplitudes of the magnetic field components are 
rescaled by the same factor. 

5. SUMMARY AND DISCUSSION 

We have argued that acceleration of AGN jets to relativistic 
velocities on scales that are much larger than the gravitational 
radius of the central black hole is most plausibly explained in 
terms of magnetic driving. This mechanism involves acceler- 
ation by the gradient of the azimuthal magnetic-field pressure 
and is distinct from centrifugal acceleration, which is often 
considered to be the dominant driving mechanism of non- 
relativistic jets. Centrifugal driving takes place in the sub- 
Alfvenic flow regime and accelerates the gas to a poloidal 
speed that is of the order of the initial Keplerian speed in the 
underlying disk. In comparison, magnetic pressure-gradient 
acceleration occurs over a much more extended region (up 
to the modified fast-magnetosonic surface) and can produce 
a much higher (relativistic) terminal speed depending on the 
initial Poynting-to-mass flux ratio /ic 2 . In the trans- Alfvenic 
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FIG. 2. — Same as Fig.0 but for the application to the superluminal jet in 3C 345. Here fDA.in = 2.7 X 10^* pc and fDA.out = 150 03a. in = 4. 1 X 10 2 pc. 



relativistic-MHD solutions presented in this paper, the termi- 
nal Lorentz factor is 7^ sa /i/2, corresponding to a rough 
equipartition between the asymptotic Poynting and kinetic- 
energy fluxes. 6 These solutions are also characterized by 
strong magnetic collimation (with the streamlines tending 
asymptotically to cylinders) and are thus consistent with the 
narrow opening angles inferred in AGN jets. 

Thermal effects could in principle contribute to the 
acceleration even in jets where a large-scale mag- 
netic field guides the flow. As discussed in VK, 
there are in general two thermal force densities: the 
pressure gradient -VP and the "temperature" force 
-l 2 po(V- VOV = -(~f 2 V/c 2 )V ■ VP. These forces ac- 
celerate the flow to 7 w . In cases where > 1, the thermal 
acceleration takes place in the nonrelativistic regime and is 
terminated by the time the speed increases to V3C s .i — i.e., 
just beyond the sonic surface that typically lies very close 
to the origin. As pointed out in § |3] this situation applies 
to proton-electron outflows; therefore, to the extent that 
AGN jets have a dynamically dominant proton component 
(as is often inferred to be the case), their acceleration to 
relativistic speeds will not be significantly influenced by 

6 Note that the Lorentz factor on the classical fast-magnetosonic surface 
is only ~ fj, 1 / 3 (e.g., Camenzind 1986), so most of the acceleration in these 
solutions occurs in the super-fast regime. 



thermal effects. 7 

Although we only considered two specific applications — 
sub-parsec-scale acceleration involving mod erate ly relativis- 
tic speeds in a radio-galaxy jet (NGC 6251; § 14. 1> andparsec- 
scale acceleration involving highly re lativ istic speeds in a 
superluminal radio quasar (3C 345; § 14.21 — there are al- 
ready several other reported cases of relativistic AGN jets 
that show evidence for a parsec-scale acceleration. In some 
cases this has been deduced from an increase in the ap- 
paren t speed of a particular superluminal component (e.g. 
Houah, Zensus, & Porcas 1996). In other cases, where there 
are observations of several superluminal components, it was 
found that the innermost one typically exhibits the smallest 
proper motion, with more distant com ponents indica ting an 
acceleration on parsec scales ( e.g..lHom an et al.| 2001l) . In the 
case of the quasar 3C 279 jet, Pin er et all (12003) inferred an 
acceleration from 7 = 8 at r < 5.8 pc to 7 = 13 at r « 17.4 pc 
using a similar approach to the one that had been employed by 

7 In the case of an electron-positron outflow, or under optically thick con- 
ditions when radiation pressure contributes strongly to the specific enthalpy, 
one can have ^> 1 . In this case the thermal driving (dominated by the "tem- 
perature" force) acts well beyond the sonic surface and accelerates the flow 
to a highly relativistic speed (see §Q. If, in addition to 3> 1, M ^ £i also 
holds, then magnetic driving takes over at the end of the thermal acceleration 
zone and eventually increases the Lorentz factor to ~ (see § [3J: this is the 
behavior obtained in the GRB jet models referenced in footnotel2l 
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lUnwin et al 1 JT997ll in 3C 345. It is also worth noting in this 
connection that a variety of observations indirectly support the 
magnetic acceleration picture for AGN jets. For example, the 
parsec-scale helical field morphology implied by our model 
is consist ent with VLBI polar ization maps of BL Lac ob- 
jects (e.g., Gabuzda et al. 2000) a nd with circular-p olarization 
measu rements of blazars (e.g., iHoman. Attridee. & Wardlel 
120011) . However, we defer a more detailed discussion of the 
additional observational implications of this model to future 
publications in this series. 

In conclusion, we reemphasize that our modeling frame- 
work is quite general and is potentially applicable to relativis- 
tic jets in a variety of astrophysical settings. In our previ- 



ous application to GRBs, the model could account for the in- 
ferred values of 7^ and of the upper limit (~ 10 14 cm) on 
the size of the acceleration region, but it could not be further 
constrained because the acceleration region in GRB sources 
is not resolved. In contrast, the motion of the radio compo- 
nents in certain microqua sar jets has been monitored on scales 
~ 10 cm (e.g. JFender"etal] ["l999). It would thus be inter- 
esting to search for evidence of extended acceleration in these 
sources, in analogy with the situation in AGN jets. 
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